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Abstract. We address the sensitivity of quantum mechanical time evolution by 
considering the time decay of the Loschmidt echo (LE) (or fidelity) for local 
perturbations of the Hamiltonian. Within a semiclassical approach we derive analytical 
expressions for the LE decay for chaotic systems for the whole range from weak to 
strong local perturbations and identify different decay regimes which complement 
those known for the case of global perturbations. For weak perturbations a Fermi- 
golden-rule (FGR) type behavior is recovered. For strong perturbations the escape- 
rate regime is reached, where the LE decays exponentially with a rate independent of 
the perturbation strength. The transition between the FGR regime and the escape- 
rate regime is non-monotonic, i.e. the rate of the exponential time-decay of the LE 
oscillates as a function of the perturbation strength. We further perform extensive 
quantum mechanical calculations of the LE based on numerical wave packet evolution 
which strongly support our semiclassical theory. Finally, we discuss in some detail 
possible experimental realizations for observing the predicted behavior of the LE. 



PACS numbers: 05.45.Mt, 03.65. Sq, 32.80.Lg 



| Present address 

§ Present and permanent address 



Loschmidt echo for local perturbations 



2 



1. Introduction 

One of the most prominent manifestations of chaos in classical physics is the 
hypersensitivity of the dynamics to perturbations in the initial conditions or 
Hamiltonian. That is, two trajectories of a chaotic system launched from two 
infinitesimally close phase-space points deviate exponentially from each other; so do 
the trajectories starting from the same point in phase space, but evolving under slightly 
different Hamiltonians. In a quantum system it is natural to consider |(0i|02)| 2 as a 
measure of "separation" of two quantum states \<px) and \4>2)- The unitarity of quantum 
propagators renders the overlap of any two states of the same system unchanged in 
the course of time. Thus, quantum systems are said to be stable with respect to 
perturbations of the initial state. However, a perturbation of the Hamiltonian can 
(and usually does) result in a nontrivial time dependence of the wave function overlap, 
suggesting a viable approach for describing instabilities and, therefore, for quantifying 
chaos in quantum systems. 

Peres [1] proposed to consider the overlap 



of the state e~ tHt ^ h \(j)o) , resulting from an initial state |0o) after evolution for a time t 
under the Hamiltonian H, with the state e~ tHt ^ h \(j)o) obtained from evolving the same 
initial state through t, but under a slightly different (perturbed) Hamiltonian H. He 
showed that the long-time behavior of 



depends on whether the underlying classical dynamics is regular or chaotic. 

In the field of quantum computing M(t) is an important concept, usually referred 
to as fidelity [2j. Moreover, M(t) can be also interpreted as the squared overlap of the 
initial state |0o) and the state obtained by first propagating |0 O ) through time t under 
the Hamiltonian H, and then through time —t under the perturbed Hamiltonian H 
(or —H from t to 2t). This time-reversal interpretation constitutes a description of the 
echo experiments that have been performed by nuclear magnetic resonance since the 
fifties [3]. When the Hamiltonian H describes some complex (many-body or chaotic) 
dynamics M(t) is referred to as Loschmidt echo (LE) [I], and this is the terminology we 
will adopt. 

By construction, the LE equals unity at t — 0, and typically decays further in 
time. Most of the analytical studies so far addressed the quantity M(t) corresponding 
to the LE averaged either over an ensemble of initial states, or over an ensemble 
of different perturbed (and/or unperturbed) Hamiltonians. M(t) has been predicted 
to follow different decay regimes in various chaotic systems with several Hamiltonian 
perturbations [HIE]. Depending on the nature and strength of the perturbation, H — H, 
one recognizes the perturbative Gaussian [USE], the non-diagonal or Fermi- golden-rule 
(FGR) [3 [71 El [TO] and the diagonal or Lyapunov [5J [H] regimes. Here, 'diagonal' and 
'non-diagonal' refer to the underlying pairing of interfering paths in a semiclassical 



0(t) = {(j) \e im ' n e~ iHt ' h \(j) ) 



(1) 



M(t) = \0(t)\ 2 



(2) 
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approach, see Sec. [2J The perturbative, FGR and Lyapunov regimes, listed above 
in the order of the (properly defined) increasing perturbation strength, constitute the 
framework for classification of LE decay regimes [6j[12]. It is important to mention that 
the full variety of system- and perturbation-dependent decay regimes is rather rich, and 
extends far beyond the above list: double-Lyapunov p2], super-exponential [2] and 
power law [15] decay regimes serve as examples. We further note that analytical results 
for the time decay of the unaveraged LE, M(t), are currently available only for very few 
chaotic systems [16J. 

The discovery [5] of the Lyapunov regime for the decay of the averaged LE in 
classically chaotic systems, M(t) ~ exp(— At) with A being the average Lyapunov 
exponent, provided a strong and appealing connection between classical and quantum 
chaos: it related a measure of instability of the quantum dynamics, such as the LE, 
to a quantity characterizing the instability of the corresponding classical dynamics, i.e. 
the Lyapunov exponent. This result awoke the interest on the LE in the quantum 
chaos community. The Lyapunov regime has been numerically observed in several two- 
dimensional chaotic systems, i.e. in the Lorentz gas [Til [T7], the Bunimovich stadium 
[T8] , the smooth stadium billiard [T9] , a Josephson flux qubit device [20], as well as in 
one-dimensional time- dependent Hamiltonian systems [TJ. 

The theory of the Lyapunov decay of the LE mainly relies on the following two 
assumptions: (i) the validity of the structural stability arguments (supported by the 
shadowing theorem [21]), and (ii) the global nature of the Hamiltonian perturbation. 
The first assumption guarantees a unique one-to-one mapping of trajectories of the 
unperturbed system to those of the perturbed system. This mapping allows for efficient 
pairing of the trajectories of the unperturbed and perturbed system in the diagonal 
approximation [22J. The second assumption implies that the Hamiltonian perturbation 
affects every trajectory of the system, and, therefore, all trajectories are responsible for 
the decay of the LE. However, this is by no means the most general situation when we 
consider possible experimental realizations of the LE. 

In the present work we extend the semiclassical theory of the LE by lifting the 
second of the two above-mentioned assumptions, i.e. we allow for a local perturbation 
in coordinate space. In this context the LE decay was previously addressed in the 
case of a strong local perturbation [23], i.e. for a billiard exposed to a local boundary 
deformation much larger than the de Broglie wavelength. Analytical and numerical 
calculations yielded a novel LE decay regime, for which M(t) ~ exp(— 2yt) with 7 being 
the probability (per unit time) of the corresponding classical particle to encounter the 
boundary deformation. 7 can also be viewed as a classical escape rate from a related 
open billiard obtained from the original (closed) one by removing the deformation- 
affected boundary segment. In this work we explore all possible strengths of a local 
perturbation and describe the transition from the weak to the strong perturbation 
regime that completes the previous picture. In particular we show that the rate of 
the exponential decay of the LE oscillates as a function of the perturbation strength as 
the escape-rate regime is approached. 
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The paper is organized as follows: In Section [2] we develop a comprehensive 
semiclassical approach of the LE decay due to local Hamiltonian perturbations of 
increasing strength. We perform a systematic analysis of the different decay regimes 
and establish their relation to the previously known decay regimes in the case of a 
global perturbation. In Section [3] we validate our semiclassical theory by numerical 
simulations. In Sec. H] we outline possible experimental realizations and focus on the 
possibility of introducing a local perturbation in microwave-cavities and ultra-cold atom- 
optics billiards. We provide concluding remarks in Sec. [5] and point to the similarities 
and differences with respect to another approach to the LE for local perturbations. 
Technical aspects of the calculations are relegated to the appendices. 



2. Semiclassical approach 



2.1. Wave-function evolution 

We address the time evolution of the wave function that describes a quantum 
particle moving inside a classically chaotic two-dimensional billiard (corresponding to a 
Hamiltonian H). We assume that initially (at time t — 0) the particle is in a coherent 
state 



exp 



h P0 



r J 



2a 2 



(3) 



Here a quantifies the extension of the Gaussian wave packet, while r and po are the 
initial mean values of the position and momentum operators, respectively. We further 
define the (rescaled [21]) de Broglie wavelength of the particle as X = h/po. 

In our description of the time evolution of the wave function we rely on the 
semiclassical approximation [25] of the wave function at a time t, 



r) = J dv' /Q(r,r',t)0 o (r') 



3(r,r',t) 



Here 



fC s (r,r',t) 



VP's 
2irih 



exp 



r,r,t) -i— 



(4) 



(5) 



is the contribution to the Van Vleck propagator associated with the classical trajectory 
s(r, r', t) leading from point r' to point r in time t. Ss{r,r',t) denotes the classical 
action integral (or the Hamilton principal function) along the path s. In a hard-wall 
billiard S$(r, r',t) = (m/2t)L?(r, r'), where L s (r, r') is the length of the trajectory s, 
and m is the mass of the particle. In Eq. ([5]), Dg = \ det(— d 2 S$/drdr')\, and the Maslov 
index v~ s equals the number of caustics along the trajectory s plus twice the number of 
particle- wall collisions (for the case of Dirichlet boundary conditions). 

Since we assume that the initial wave packet is localized around r within a, only 
trajectories starting at points r' close to r are relevant for our semiclassical description. 
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Figure 1. Sketch of a typical trajectory s (full line), connecting a point r' (within 
the circular extension of radius a of the initial wave packet) to the point r, where the 
evolved wave function is evaluated, together with the central trajectory s (dashed line) 
that reaches the same final point, but starts at the center r of the wave packet (with 
momentum p s ). The linearization of Eq. ((6]), together with the conditions discussed in 
the text, allow to represent all the trajectories s contributing to Eq. (j4|) by the single 
reference trajectory s. 



Thus we can expand the action integral Ss(r,r',t) in a power series in (r' — r ). In 



Appendix A we show that the power series can be terminated at the linear term, 

S»(r, r', t) « S s (r, r , t) - p s ■ (r - r ) , 
if the wave packet is narrow enough, so that the condition 



a < 



h 



(6) 



(7) 



1/X + l/o- 

is satisfied, where Zl = Po/ m ^ is the Lyapunov length. In Eq. (El), s(r, r ,t) is the 
central reference trajectory into which s(r, r', t) gets uniformly deformed as r' — > r , and 
p s = —dS s (r,r ,t)/dr denotes the initial momentum of the trajectory s (see Fig. CQ. 



Substituting Eq. ([6]) into Eq. 
obtain 



J5]), and performing the integration in Eq. (j3J) we 



Mr 



2nh 



E 

s(r,ro,t) 



/C s (r,r ,t) $ (Ps 



with 



$o(p) = 



dr 

2nh 



exp 



^o(rj 



(8) 



(9) 



a 



exp 



a 



2h 2 



,P - Poj 



the momentum representation of the initial wave packet. 

We now consider a related billiard, corresponding to the perturbed Hamiltonian H, 
that differs from the original (unperturbed) billiard by a deformation of the boundary 
segment B\ of width w (see Fig. [2]). The perturbation is thus local, and will be 
characterized by its extent (depending on the ratio between w and the cavity perimeter 
P) and its strength (that will be quantified in the sequel). In view of Eq. (JHD, the wave 
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Figure 2. Sketch of our model system of a particle moving inside a chaotic billiard 
of area A. The perturbation consists of a deformation localized in a region B\ (of 
width w) of the billiard boundary. The complementary set Bq of the boundary is 
unaffected by the perturbation. The two trajectories s (red solid line) and s (red 
dashed line) starting from ro with different momenta correspond, respectively, to 
the unperturbed and perturbed Hamiltonian. The diagonal approximation entering 
Eq. ([12]) identifies both of them and assigns an action difference given by Eq. (I14| . 
The starting momentum of the solid red trajectory belongs to the set V\. The third 
trajectory (blue solid line) hits the boundary only at Bo and therefore is the same 
for both the unperturbed and perturbed systems. Hence the action difference of the 
corresponding trajectory pair is zero. The starting momentum of the blue trajectory 
belongs to the set Vq. 



function describing the evolution of the particle (starting from the same initial state </>o) 
can be written as 

Mr) = 2nh J2 /C s ~(r,r ,t)$ (p s ~). (10) 

s(r,ro,t) 

The sum now runs over all possible trajectories s(r, ro, t) of a classical particle that 
travels from tq to r in time t while bouncing off the boundary of the perturbed billiard. 



2.2. Wave-function overlap for local perturbations 

According to Eqs. (!H|)- (|T0|) and the definition (pQ) of the LE amplitude, we have 



0(t) 



dr 



(11) 



A 



J A r. 



S S 
% 



x exp 



%{Vs)%{Ps) 



where A stands for the billiard area. The shadowing theorem [21] allows us to employ 
the diagonal approximation (s ~ s) in the case of a classically small perturbation 
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thus reducing Eq. (ITTj) to 



0{t) 



dr 



A 



E 

?(r,r ,i) 



_D S exp 



-AS a (r,r ,t) 



Wo (P. 



where 



W (p) = |$ (p) 



^2 



exp 



PoJ 



(12) 



(13) 



is the probability distribution of the particle momentum. In billiards the action 
difference between the two trajectories traveling between the same initial and final points 
in the same time t can be written, in terms of their length difference AL S , as 



AS s (r,r ,t)=S s -S s 



(14) 



2m 



Ps 



Ps ~ Ps 



t = p s AL s (r,r ,t) 



2m m 

Using the Jacobian property of the Van Vleck determinant D s , we can replace the 
integral over final coordinates in Eq. (fl2"l) by an integral over the initial momenta and 
obtain 



0(t) = / rfp^ (p)exp 



h 



p AL(r ,pt) 



(15) 



The dependence of AL on the product pt stems from the fact that in billiards, changing 
the magnitude of the momentum only modifies the traveling time, but does not affect 
the path. 

We now introduce a sequence of momentum sets P„(r , t), with n = 0, 1, 2, . . ., such 
that for any p € V n the classical trajectory, starting from the phase-space point (r , p), 
arrives at a coordinate-space point r G A after time t while visiting the deformation- 
affected boundary segment (£>i in Fig. [2]) exactly n times. Thus, the trajectories with 
the initial momentum in Po( r o, t) only undergo collisions with the part of the boundary 
unaffected by the deformation (£> in Fig. [2]) rendering AL = 0, so that 



0(t)=J20n(t) 



n=0 



with 



O (t) 
O n (t) 



dp W (p) 



[ dp W (p) e ipAL{ro > pt)/h for n > 1 . 



(16) 



(17) 



Since we are studying billiards, the integrations over momenta can be simplified by 
working in polar coordinates (p, 6) and considering the sets Q n (r ,pt) of angles 6 such 
that p = (p, 6) G V n (r Q ,t) iff 9 G n . For a classically chaotic dynamics the set G>o 
shrinks with increasing time t, and Bo becomes the fractal set defining the repeller of 
the corresponding open (scattering) problem in the limit t — > oo. Eqs. (1171) and (|T8l) 
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can, respectively, be written as 



o (t) = dpp de w„(p) , (19) 

Jo Joq 

poo /» 

o n (t) = dpp de w (p) eip^o^m/h for n > 1 _ ^ 20 ) 

</o </e„ 

For long times t, where many trajectories contribute to the semiclassical expansions, 
the angular integrals over n can be replaced by integrals over all angles, weighted with 
the measures p n of the corresponding sets. Assuming ergodicity one readily obtains 
an approximation for the probability for a trajectory of length I = pt/m to visit the 
boundary region B\ exactly n times: 

"^"^{rS^H)' (21) 

where Id is the average dwell length of paths in the related open chaotic billiard obtained 
from the original (closed) one by removing the boundary region B\. This corresponds 
to the classical escape rate of the open cavity that for particles with momentum po is 
given by 

= < 22 > 

For a chaotic cavity with an opening w much smaller than its perimeter P we can 
approximate [27] Zd ~ ttA/w, and therefore 

l^~ A - (23) 

In our case the escape rate 7 yields a measure of the perturbation extent. The classical 
escape rate of an open cavity controls the fluctuations of the transmission coefficients, 
and therefore approximations such as ( 1231) have been thoroughly examined in the context 
of quantum transport [28J. 

According to the previous discussion we can approximate O n (t), with 71 = 0,1,2,..., 
by the averages 

O (t)= dpp Po (pt/m) d9W (p), (24) 
Jo Jo 

poo p'2n 

O n (t) = / dpp(e ipAL ^ pt ^ h ) e P n(pt/m) / d9 W (p) ,n>l. (25) 
Jo ™ </o 

The mean value ( . . . )e n should be taken over the set n (r o , pt). The chaotic nature of 
the dynamics will enable us to treat the averages over G n in a statistical way. In view 
of Eq. (QJ), the ^-integral in Eqs. flMD and ([25]) yields 

r 27T 2cr 2 T a 2 1 /2o- 2 \ 

d6Wo(p) = -^exp --(p 2 +pl) I ^— Po p\, (26) 

where Jo denotes the modified Bessel function. 

As usually assumed in the Loschmidt echo studies, we restrict our analysis to 
"semiclassical" initial wave packets <po(r) with sizes much larger than the de Broglie 
wave length, 

A < a . (27) 
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This assumption, together with condition (j7j), defines the interval for the dispersion 
a, where the semiclassical approach is reliable, and hence yields restrictions to the 
parameters of the billiard. Employing condition (12?]) enables us to use the asymptotic 
form Io(x) ~ e x /\ / 2nx, valid for large x, in Eq. fl26l) . Thus, the probability distribution 
function for the magnitude of the initial momentum is given by 



p d0W o (p) « \J^- exp 
Jo n V *Po 



'2 



(2f 



We note that Eq. (j2"Bj) provides a good approximation to the exact distribution function 
already for a > 2A. Under this assumption the p-integrals in Eqs. ( |24|) and ( 1251) are 
dominated by the contributions around po, and we can write (in view of Eq. (l2Tj) ) 

Oo^-e" 7 *, (29) 
OJt) « M! e -7* / e i P oAi(r , Po t,e)/a\ for n > x _ (3Q) 

Since in Eq. (|30|) all classical quantities are evaluated for an initial momentum with 
magnitude po, the mean values ( ... }e n should be taken over the sets Q n (r Q ,p Q t). 
However, for long times and a chaotic dynamics we do not expect these mean values 
to depend on r . In the next section we will further invoke the chaotic nature of the 
underlying classical dynamics in order to estimate the mean values and therefore the 
LE average amplitude. 



2.3. Averages over trajectory distributions 

For classically small perturbations the action, respectively, length difference (Eq. (TUft ) 
between a trajectory s (solid red segment in Fig. [2]) and its perturbed partner s (dashed 
segment) is given only by the contributions accumulated along the n encounters with B\. 
Differences in length arising from the free flights between collisions with the boundary 
(Bo + Bi) are of higher order in the perturbation strength and will not be considered. 
We can then write 

■n 

AL = £; (31) 

3=1 

where the deformation function u("&j,£j) is the length difference accumulated in the 
j-th collision with B±, depending on the impinging angle i?j G (— 7r/2,7r/2) and on the 
coordinate £j G (0, w) of the hitting point. The number n of collisions with B\ is a 
fraction of the total number of collisions pot/mlf. The mean bouncing length if can be 
approximated by nA/P [29], and we suppose U 3> k since P ^> w. We note that for 
small perturbations AL depends on p only through n. 

Given the chaotic nature of the classical dynamics and the fact that the collisions 
with B\ are typically separated by many collisions with Bq we assume {$j} and 

to be random variables. Furthermore, assuming a perfect randomization of the 
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trajectories within the billiard the probability distribution functions for and 
are, respectively, 

PM = ^r and p d0 = ~- (32) 
Then, treating the random variables as uncorrelated we obtain 



e 



ipoAL(ro,pot,9)/h\ _ /„iu/X 



>e„ = ¥"'T • (33) 



'Or, 

where the average (...) is defined as 

/7r/2 pw 
aimm / dtPtwmt) (34) 
-tt/2 J0 

for a function f Once we specify the shape of the perturbation, the average in 

the right-hand side of Eq. (|33|) is readily calculated from the probability distributions 
of Eq. (|32l) . For instance, for a piston-like deformation (see Appendix B ), 

u(#,f) = 2/icos#, (35) 

so that the average reads 

(e m/x ) = 1 - | [H X (2VA) - iJi(2ty*)] , (36) 

where Hi stands for the Struve H-function of the first order, and J± is the first order 
Bessel function of the first kind. Here we note that for the case of the piston-like 
boundary deformation the ratio h/X serves as a measure of the perturbation strength, 
whereas w (and therefore 7) quantifies the extent of the perturbation. 

At this stage, however, we will keep our discussion general and do not specify the 
details of the local perturbation. The substitution of Eq. (133]) into Eq. (1301 yields 
the average LE amplitude. The latter is usually not an observable quantity, however 
it will be helpful towards our semiclassical calculation of the LE. The condition of a 
classically small perturbation that we have adopted throughout our work, implies that 
y/ (u 2 ) <C w P (which for the case of the piston-like deformation is equivalent to 
h w P). Quantum mechanically the perturbation can be characterized by the 
extent w and a deformation strength 

X = ^Z. (37) 

For x C 1 we will be in the quantum perturbative regime [3, El E] , which will not be 
considered in this work. Increasing the deformation strength x we anticipate a richer 
variety of regimes than for the case of the LE under global perturbations [5J [6] since the 
perturbation extent enters as another relevant parameter. 



2.4- Loschmidt echo for local perturbations 

According to Eqs. (j5J) and ffTTl) the semiclassical expansion for the LE contains terms 
involving four trajectories. The diagonal approximation, leading to Eq. ( 1T21) for the 
LE amplitude, reduces the LE to a sum over pairs of trajectories. Consequently, the 
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semiclassical form of the LE must take into account the different possibilities for each 
trajectory of the pair to hit (or not) the region of the boundary where the perturbation 
acts. We can therefore decompose the LE as 

M(t) = M nd {t) + M d {t) , (38) 

where we have introduced the non-diagonal and diagonal contributions according to 

M nd (t) = J dpW (p)e- tpAL{ro > pt)/h J_ dp'W (p')e ip ' AL{ro > p ' t)/h (39) 

and 

M d (t)= JdpW (p)£ dp'Wo(p') (40) 

x exp |^ \p'AL{r , pt) - P AL(r , pt)] j . 

Here, the set e p of momenta p' is defined such that two trajectories starting from the 
phase space points (r , p) and (r ,p'), stay "close" to each other in phase space during 
time t, and thus are "correlated" with respect to the perturbation; e p is the momentum 
set complementary to e p . We give a quantitative definition of e p below. Following 
the standard notation introduced in Ref. [5], we call diagonal term the one resulting 
from the identification of pairs of trajectories where the effect of the perturbation is 
correlated, that is, when p' 6 e p . In the non-diagonal term we consider the pairs of 
trajectories uncorrelated with respect to the perturbation, including the case where one 
or both orbits are unperturbed. As noted at the beginning of this section, each of the 
trajectories of the above pair already incorporates a diagonal approximation between a 
perturbed and an unperturbed trajectories with the same extreme points. 

2.5. Non-diagonal contribution to the Loschmidt echo 

Calculating the LE as an average over trajectory distributions forces us to take into 
account pairs of trajectories and the possible correlations among them. The correlations 
are particularly important for M d (t), as we show in Sec. 12. 61 On the other hand, in 
our discussion of the last chapter we established that for the calculation of M nd (t) the 
two trajectories of the pair can be considered to be uncorrelated with respect to the 
perturbation, and the averages can be performed independently. Assuming in addition 
that the measure of the momentum set e p is small compared with that of the momentum 
set effectively represented by the distribution Wo(p) we write 

2 



M nd (t) 



n=0 



(41) 

Substituting Eqs. (I29H30I) and (133]) into Eq. f|4T]) and performing straightforward 
algebraic operations we arrive at our central result for the non-diagonal contribution 
to the LE: 

M nd {t) « e~ K7i , (42) 
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where 

K = 2(l-^(e iu/x )) . (43) 

Thus, for the case of the piston-like boundary deformation the average phase factor due 
to a single visit of the perturbation region B\ by the particle is given by Eq. ([3*61) . so 
that Eq. fl43]) reads 

k = 7rHi(2/i/A) • (44) 

In subsection 12.71 we study the emergence of the Fermi-golden-rule and escape- 
rate regimes of the decay of M nd (as well as the transition between the two regimes) 
depending on the strength of the perturbation. In sections [3] and H] we discuss their 
implications for numerical simulations and possible experimental observations. 



2.6. Diagonal contribution to the Loschmidt echo 

To proceed with the calculation of the diagonal contribution to the LE, Eq. (JIT)]) , we first 
need to specify the set e p such that two trajectories of time t, starting from the phase 
space points (r , p) and (r , p' e e p ), stay "correlated" during time t. As in Sec. \'2.2\ it is 
convenient to work with polar coordinates, p = (p, 9) and p' = (p', 9'), in which the set 
e p can be defined as follows: for every p' 6 e p one has \p' — p\ < Ap and \8' — 9\ < AO. In 
turn, Ap and AO are subject to the requirement that the two trajectories stay "close" to 
each other in phase space. Indeed, any two "correlated" trajectories must have the same 
number of collisions with the billiard boundary. This condition leads to \p' — p\t/m < /f. 
Moreover, the two trajectories must also have the same number of collisions with B\. The 
spatial separation between the two trajectories at the first collision with the boundary 
is approximately given by \9' — 9\l{. The condition that after time t this separation is 
smaller that the size w of B\ is \9' — 8\l{ exp(At) < w. Here we have used the property 
that for chaotic dynamics two arbitrary, initially close trajectories deviate exponentially 
from each over with a rate given by the average Lyapunov exponent A. Thus, we can 
estimate the measure of the e p -set to be 

Ap = — - and AO = — exp(-At) . (45) 
t If 

Using this quantitative description of e p for the evaluation of Eq. fj40|) we obtain for the 
diagonal contribution to the LE 

/rp+Ap r6+A6 
d P W (p) / dp'p' / d9'Wo(p') (46) 
Jp~Ap Je-Ae 

x exp | l - \p'AL(r , p% 6') - P AL(r , pt, 9)] 

We now argue that for boundary deformations of moderate strength the exponent in 
the integrand on the right hand side of Eq. (1461) can be neglected. The argument of the 
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exponent is given by the total differential of the function pAL(r Q ,pt, 
its absolute value can be bounded by 

(p'-p)AL p(6' -6)dAL 



h h 89 

Here we have used that, as discussed in Sec. 

yt) of collisions with B\. Then 



< 



ApAL 



h 



pA6 dAL 



h 89 



and therefore 



(47) 



number n 



ApAL 



h 



n id 



AL is independent of p for a fixed 



(4f 



and 



pA6 dAL _ pA6 



du ddj 
Mj d6 



+ 



du 8£j 



A6 Ce 



xt 



w 



C 



(49) 



r = -l d — 

XXd-d 



(50) 



where we have introduced the dimensionless quantity 

/f I du 
"AW, 

This implies that the exponent in Eq. (j4"6l) is smaller than unity if x ^ ^d/k an d 
C < lf/w. Since the ratios l^/k and lf/w are assumed to be large, the above inequalities 
hold for a wide range of deformations. (Note that for the piston-like deformation, 
Eq. f[3"5j) . x — (8/3) 1 / 2 h/X and C = 0, and hence the exponent is small if h < Xl^/k-) 
Neglecting the exponent in Eq. fj4*6|) we obtain 



M d (t) 



p+Ap p0+A0 

dpW (p) I dp'p' / d0'W o (p') 

p-Ap J6-A6 

AApAO I dppW 2 (p) . 



(51) 



In the second line of this equation we have taken into account that Ap p and 
AO 1 for times t much longer than the dwell time t^- Under the assumption fT27|) of a 
"semiclassical" initial wave packet, the integral over p in Eq. (J5T1) is dominated by the 
contribution around p and we find 

M d (t)^^mY- (52) 



We note that the exponential dependence of M d on At does not contain the 
perturbation and arises from a classical probability distribution, as in the standard 
Lyapunov regime [5]. The dependence with respect to the perturbation extent, w, 
appears in the prefactor. 



2. 7. Decay regimes of the Loschmidt echo 

According to Eq. (1381) the full LE M(t) is the sum of the non-diagonal and diagonal 
contributions, M nd (t) and M d (t), given by Eqs. (j4"2"|) and (|52|) respectively. We first 
argue that, unlike in the case of a global Hamiltonian perturbation, the non- diagonal 
contribution will typically dominate over the diagonal term. The most favorable regime 
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to observe the diagonal term would be that of perturbations satisfying A < ivy. 
Therefore, the necessary requirement for resolving the diagonal term is A < 4y. Since 
A — 7 = /iks > (with Hks the Kolmogorov-Sinai entropy of the chaotic repeller [30J) we 
see that M d (t) can possibly prevail over M nd (t) only if /iks < 3A/4. Taking into account 
that /iks = A for a closed system, we see that only relatively open cavities would allow 
to observe M d (t). Hence for further analysis of the LE decay we mainly focus on the 
non-diagonal contribution M nd (t). 

According to Eq. ( T42|) the LE displays different decay regimes depending on the 
strength of the perturbation. Thus, for weak perturbations characterized by x ^ 1 one 
can expand the phase factor e %u l x in the Taylor series to the second order in u/A to 
obtain 9ft(e m / A ) « 1 — % 2 /2 and, therefore, 

M nd (t)^e~ xV for X <1 (FGR). (53) 

The rate of the exponential decay given by Eq. (J53"|) depends on the perturbation 
strength \i m analogy to the Fermi-golden-rule regime found for global perturbations, 
but is dressed with 7 that provides a measure of the fraction of phase-space affected 
by the boundary deformation. On the other hand, in the limit of strong perturbations, 
X 1, the LE decay rate is independent of the perturbation strength \ an d is entirely 
determined by the extent of the deformation quantified by 7. Indeed, (e m / A ) — » as 
X — > 00 leading to 

M nd (t) « e" 27 * for x > 1 (escape-rate) . (54) 

This is the escape-rate dominated decay regime previously reported in Ref. [23]. We 
finally emphasize that Eq. (j4"2"l) . and therefore Eqs. (153]) and fl54"j) . hold only if the 
conditions ([7]) and ( 1271) are satisfied, i.e. 

A < a < a/^l , (55) 

where /l is the Lyapunov length. 

As a particularly interesting feature, the transition between the FGR and escape- 
rate regime is non-monotonic, i.e. the LE decay rate Ky, in general, oscillates as a 
function of the perturbation strength x while approaching the asymptotic value 27. 
Figure [3] illustrates these distinct oscillations for the case of a chaotic billiard with 
the Hamiltonian perturbation generated by a piston-like boundary deformation, see 



Appendix B In other words, for a fixed time t and starting in a minimum, the LE 



can increase by orders of magnitude upon varying the perturbation strength x- The 
strength of the oscillations depends on the geometry of the boundary perturbation. The 
oscillations are particularly pronounced for the piston-type geometry where they only 
very slowly merge into the escape-rate limit. 
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Figure 3. Decay rate (in units of 7) of the Loschmidt echo in a chaotic billiard 
(Eq. (j43|) ) as a function of the perturbation strength x- We have chosen a piston-like 
boundary deformation (see | Appendix B| for which the LE decay follows Eq. (|44ll and 
X = (8/3) 1 / 2 h/X. The Fermi-golden-rule decay regime (Eq. (|53|) ) is recovered for x ^ 1 
(red dashed line). As x ~ * 00 the decay rate asymptotically approaches the value 2 
(black dashed line) representing the escape-rate regime (Eq. ([54]) ). 

3. Numerics versus semiclassical predictions 

3.1. Finite-size corrections to the semiclassical limit 

In deriving our approximate expression for the non-diagonal contribution to the LE 
(Eq. (H2l) ) we made the following two important assumptions. Firstly, we restricted 
our discussion only to those initial wave packets whose dispersion a is much larger 
than the de Broglie wavelength X, see Eq. (|27|) . In numerical simulations, however, one 
can only address finite cr/X ratios, so that the theory must be improved accordingly 
to be capable of accounting for the results of the simulations. Secondly, in Sec. [2] we 
used the simple expression, given by Eq. (121]) . for the probabilities p n (l) for a classical 
trajectory of length / to visit the deformation region n times. While this expression 
well approximates p n (l) for I > Ij an improved formula is needed to correctly treat the 
numerically assessable range of lengths, I ~ U- In this subsection we address these two 
important issues and present an expression for M nd (t) appropriate for a quantitative 
comparison with the results of the numerical simulations. 

The non-diagonal contribution to the LE, as given by Eq. ( 14T1) . is given by the 
square of the absolute value of the sum of overlaps O n {t) (with n = 0,l,2,...) defined 
according to Eqs. ( 124"1) and (1251) . Now instead of using the asymptotic form of the Bessel 
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function I (like in Sec. \2.2\i we keep the analysis general by writing O n (t) as an integral 
over the dimensionless momentum variable z = p/po'- 



and 



O (t) 



O n (t) 



2a 2 



2a 2 
A 1 " 



dz z po(zpot I m) exp 



a 



2a 2 



dzz (e izu/x ) n 



(56) 



(57) 



x Pn{zpot/m) exp 



O" 2 , 9 



2a 2 



n > 1 



We now address the probability distributions p n (l). The central building block 
of our analysis here is the probability g(l)dl for a classical trajectory with the length 
between (L + /) and (L + I + dl) to end on the boundary deformation region B\ subject 
to the condition that the previous encounter with E\ took place at length L. Assuming 
ergodicity for the billiard system under consideration we approximate g by the Heaviside 
step function 9 as 



0(0 »f 0(1 -Jo), 



(58) 



where the length Zo is the (average) minimal length that a trajectory starting on B\ must 
have to return to B%. In view of Eq. (1581) one readily obtains the following approximation 
for the survival probability: 

p (l) « 9(l - Z) + 9(1 - Z ) e-C-WA" . (59) 



Then, the visit probabilities for n > 1 are calculated as 



Pn(0 



dZ n _i . . . / dhp (l - Z n ) 
./o 

x #(Z„ - Z ft _i) p (^n - Z„-i) . . . #(Z 2 - li) Po(^2 - Zi) g{k) Po(h) 
and can be shown to satisfy the following recursion relation for n > 1, 



(60) 



A 



,(0 = 9(1 - nlo) 9((n + l)/ -1)1-^ p fc (Z) 



(61) 



fc=0 



+ 0(/-(n+l)Z„) 



n-l 



J^p fc ((n + 1)Z 



k=0 



n-l 
fc=0 



p k ((k + l%) n-(n + l% 
(n-k)\ V 



n—k 



exp 



Z-(n + l)Z c 
Zh 



Note that Eq. (IBTI) together with f|59l) simplifies to Eq. (l2Tj) if one puts Zo = 0. However, 
for trajectories of lengths comparable to Id the minimal return length Zo can not be 
neglected and Eqs. (159]) and (J6"T|) must be used in Eqs. fl56l) and (j5"7|) to yield the LE. 
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3.2. Numerical simulations 

In order to support our semiclassical predictions we present in this section numerical 
quantum mechanical calculations for a local perturbation. We use the Trotter- Suzuki 
algorithm J3TJ [32] to simulate the dynamics of a Gaussian wave packet inside a 
desymmetrized diamond billiard (DDB). The DDB is defined as a fundamental domain 
of the area confined by four intersecting disks of radius R centered at the vertices of 
a square. We denote the length of the largest straight segment of the DDB by L (see 
Fig. HJ). As proved in Ref. [33], the DDB is fully chaotic and thus has been previously 
considered for studying aspects of quantum chaos [31]. 




(a) 




Figure 4. Desymmetrized diamond billiard: (a) the fundamental domain of the four- 
disk billiard, (b) the initial wave packet (with momentum enclosing an angle a with 
the horizontal) in the case of a local piston-like boundary deformation (defined by a 
width w and displacement h). 



Our semiclassical analysis is valid for an arbitrarily shaped local perturbation acting 
on a region B\ (of width w) of the boundary. A perturbation with the shape of a circular 
segment was used in Ref. [23]. In our present numerical simulations we chose a piston- 
like perturbation (Fig. 0b), for which analytical results can be readily obtained (see 
Eqs. (J35J) and (|H]L and |Appendix Bp . 

In Fig. [5] we present the LE decay for a DDB with L = 1000, R = 1311 (in units of 
the lattice spacing of the underlying tight-binding model) and a piston-like perturbation. 
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Figure 5. The Loschmidt echo decay in the DDB for tree different values of the 
piston-like deformation strength (see Fig. [4]): (i) x ~ 1-63 (h = 3 and A = 3), (ii) 
X « 1.75 (h = 3 and A = 2.8), and (hi) x ~ 5.44 (h = 10 and A = 3). The solid-line 
curves are obtained as the result of numerical simulations, whereas the dashed-line 
curves show the semiclassical predictions (see text). The initial wave packet dispersion 
is a = 8, and the other system parameters are L = 1000, R — 1311 and w = 120. The 
time is given in units of the dwell time 1/7, with 7 defined in Eq. (|23p . 

For the present geometry the dwell length is Id = (P/w)l{ ~ 18.7 If, with If being the 
mean free flight path. The initial wave packet has the dispersion o = 8; its momentum 
direction is chosen to be parallel to the longest straight segment of the DDB (a = 
in Fig. Hlb), but we have verified that the LE decay rate is independent of a [35]. The 
three numerical (solid-line) curves in Fig. [5] correspond to the following values of the 
deformation strength: (i) x ~ 1-63 (h — 3 and X = 3), (ii) x ~ 1-75 (h — 3 and X = 2.8), 
and (iii) x ~ 5.44 (h = 10 and X = 3). These numerically obtained LE curves decay 
almost exponentially for times up to '-ft w 4 before turning over to a regime with strong 
irregular fluctuations around a saturation value pp. The three curves shown illustrate 
the non-monotonous dependence of the decay rate with the perturbation strength. 

The corresponding semiclassical LE decay curves - the three dashed-line curves in 
Fig. [5] - were obtained by doing the integrals in the right-hand side of Eqs. ( 1561) and 
( 1571) numerically, with p n (l) probability distributions determined in accordance with 
Eqs. (|59~1) and (|6T|) ; the minimal return length was taken to be Iq = 0.16 1&. The infinite 
sum in the right-hand side of Eq. (jUJ) was terminated at n = 8. 

The good agreement between the semiclassical predictions and the results of the 
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full quantum-mechanical computation is evident in Fig. The fact that the obtained 
LE decay deviates from the purely exponential one is entirely due to the fmiteness 
of the a/X ratio: the analytical results of Sec. [2] are recovered in the limit given by 
Eq. (I55I) . which proves challenging in numerical simulations [36] . As we discuss in the 
next section, however, this limit is naturally achievable in laboratory experiments with 
ultra-cold atom-optics billiards, so that the latter provide a viable model system for 
studying the LE from local Hamiltonian perturbations. 

4. Experimental realizations of Loschmidt echo with a local perturbation 

Experiments on the LE are of foremost importance since they render crucial information 
about quantum dynamics of physical systems and their decoherence mechanisms [37J. 
While the examples discussed in the introduction show that the agreement between 
the semiclassical theory of the LE and numerical simulations is quite successful, the 
situation is less satisfactory concerning experiments. 

LE experiments were first performed on nuclear spins of organic molecules using 
NMR techniques [31 [38] - The decay of the polarization was found to be quite insensitive 
to the coupling to external degrees of freedom or the precision of the reversal. The 
Gaussian decay of the experimentally measured LE is at odds with the one-body 
semiclassical theory, and many-body aspects of the problem have been pointed to be at 
the origin of such a behavior jTUJ, [39J EQ]. 

In Ref. [23] a principle experimental scheme for measuring the LE for local boundary 
perturbations was proposed based on a ballistic electron cavity with a small ferromagnet 
attached acting as the local perturbation. Such a setting provides a link between 
spin relaxation in a mesoscopic cavity and LE decay. Here we discuss two further 
experimental settings which appear suitable for a measurement of the echo decay: 
microwave and cold atom cavities. 

Microwave experiments allow the independent measurement of individual scattering 
matrix elements for the unperturbed and perturbed systems |41j. The cross-correlation 
of these matrix elements can then be calculated, and going into the time domain, 
the scattering fidelity is obtained. The latter is a good representation of the usual 
average fidelity amplitude when appropriate ensemble and/or energy averages are taken. 
Correspondingly, the LE can also be constructed from measurements of scattering matrix 
elements. The observation of the Lyapunov decay regime for a global perturbation 
in the microwave cavity can then be envisioned. However, reaching the required 
long time scales remains as an experimental challenge. On the other hand, the 
corresponding (escape-rate) decay regime for local perturbations might be easier to 
reach experimentally. Moreover, microwave billiards appear to be rather suited to 
investigate local boundary perturbations, since the piston-type deformation presented 
in our work can be directly realized in a microwave billiard setup. The width of the 
piston determines the exponent of the LE time decay. Furthermore, by moving the 
piston the perturbation strength \ can De directly controlled and tuned. Hence, by 
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devising sufficiently large microwave cavities to approach the semiclassical limit it seems 
promising to experimentally reach both, the escape-rate regime for large \ and the non- 
monotonic dependence of the decay rate of the LE on the perturbation strength (see 

Fig. ED. 

Studying quantum chaos in the laboratory by recreating a delta-kicked harmonic 
oscillator in an ion trap was proposed a decade ago by Gardiner and collaborators [42J, 
and a number of fruitful approaches have since then been developed and successfully 
realized using ultra-cold atoms confined to optical billiards [15] . For instance, the decay 
of quantum correlations has been measured by echo spectroscopy on ultra-cold atoms 
using the detuning of the trapping laser as a perturbation [2]. Below we focus on the 
time evolution of clouds of ultra-cold atoms in optical billiards, and show that they 
provide a viable system for experimental investigation of different perturbations and 
various regimes of the LE decay. The perturbations can be global, such as in the cases 
previously studied, but also local. Since the large-scale separation of system parameters, 
given by the condition (1551) . is attainable in these experimental systems we expect that 
a direct support for of the theoretical predictions of Sec. [2] can be obtained. 




Figure 6. Atom-optics billiard: a sketch of the focal plane of the hollow laser beam. 

In a typical microwave echo (or Ramsey) spectroscopy experiment [44J a cloud of 
ultra-cold Rb atoms is loaded into an off-resonance optical trap. For the purpose of 
our study the role of the trap can be played by a hollow laser beam with the cross 
section corresponding to the geometry of a chaotic billiard of interest. The fabrication 
of such hollow laser beams, as well as the manipulation of atoms inside them, can now be 
performed with a high level of precision [15]. The atomic cloud, after being positioned 
inside the hollow beam in its focal plane and accelerated (or "kicked") as a whole to 
a nonzero average momentum, is let to evolve freely in an effectively two-dimensional 
billiard, see Fig. El 

The Rb atoms used in echo experiments jH] are initially prepared in a quantum 
state \l/o equal to a direct product of an internal atomic state \s) and a spatial state 
described by a wave function 0o( r ), i-e. = \s) <S> 0o( r )- The internal state evolves 
in a coherent superposition of the two hyperfine sub-states, denoted by | j) and 1 1), of 
the ground state of rubidium. The | |)-component of the total wave function of an atom 
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experiences a laser field potential Vj_(r) which is, in general, different from the potential 
V|(r) exerted by the same laser on the | |)-component. The relative difference between 
the two optical potentials is given by the ratio co»hf/Al, where /ic<;hf is the energy of the 
hyperfine splitting of the ground state, and A L is the laser detuning from the frequency 
of the transition between the ground state and the first allowed excited state of the Rb 
atom. The application of a sequence of 7r/2 microwave pulses during the time evolution 
of the atoms, followed by a measurement of the populations of the | j)- and | ^-sub- 
states at the end of the evolution, allows one to determine the LE (corresponding to 
the spatial wave function 0o( r )) due to the difference between the potentials Vj_(r) and 
Vj-(r) as a function of the evolution time. 

In order to measure the LE decay due to local perturbations two different lasers 
have to be used. The first laser is to produce the confining hollow beam with the cross 
section of a desired (billiard) geometry, and has to be tuned as to exert approximately 
the same potential Vbm(r) on the both | j)- and | |)-sub-states. The beam of the second 
laser plays a role of the local Hamiltonian perturbation. It has to be placed inside 
(and aligned with) the hollow beam of the first confining laser, and its width should be 
much smaller than the linear scale of the billiard (see Fig. [6]). The frequency of this 
second laser (and perhaps its position inside the billiard) determines the perturbation 
strength x- Altering this frequency changes the difference between the potentials 5Vj_(r) 
and SV[(t), produced by the second laser and acting differently on the | J.)- and | f )- 
substates, respectively. Thus, an echo spectroscopy experiment performed in such a 
system would measure the LE decay due to the difference of the atomic potentials 
V|(r) = Vbui(r) + 8Vi(r) and Vj-(r) = Vbm(r) + SV^(r); this difference is localized in an 
area much smaller than that of the billiard. 

To date one is typically able to experimentally prepare and manipulate clouds of 
Rb atoms as cold as 1 fiK. This temperature corresponds to the thermal speed of about 
1.3 cm/s. At the same time, by first placing the atoms inside a far-off-resonance dipole 
trap then moving the trap and finally switching it off one can accelerate the atomic 
cloud as a whole up to 10 cm/s. Such a momentum kick can nowadays be easily realized 
in a laboratory, and does not significantly increase the temperature of the atoms. As 
a result one obtains a cloud of atoms moving as a whole with an average momentum 
that corresponds to a rescaled de Broglie wave length A ~ 10 ran. The dispersion of the 
atomic cloud can be shrunk to a ~ 1 /im. Under this conditions the number of Rb atoms 
composing the cloud can reach 10 5 that is well sufficient for the Ramsey-spectroscopy- 
type measurements. The hollow laser beam, producing the billiard confinement, can 
reach L « 1 cm in linear size. Assuming the Lyapunov scale Zl of the billiard to be 
of the same order of magnitude as L one arrives at the following estimate for the scale 
separation of the system parameters: 



which well satisfies the restriction given by Eq. fl55l) . Under the conditions specified 
above one can control the atoms for up to 5 s before the cloud gets significantly elongated 




(62) 
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in the axial direction of the hollow laser beam. Such a time allows a single atom to 
experience about 50 bounces with the billiard boundary, which according to our theory 
is sufficient for observing the LE decay regimes predicted in this work. 

The above considerations show that atom-optics billiards constitute promising 
candidates for experimentally investigating different regimes of the LE decay due to 
local perturbations of the Hamiltonian, and in particular put in evidence the escape- 
rate regime and the predicted oscillations of the LE as a function of the perturbation 
strength. 

5. Conclusions 

In this work we have studied the time decay of the Loschmidt echo in quantum systems 
that are chaotic in the classical limit, due to Hamiltonian perturbations localized in 
coordinate space. We have provided the corresponding semiclassical theory of the LE 
for coherent initial states evolving in two-dimensional chaotic billiards. 

In addition to the FGR decay regime, which is well-known for the case of 
global Hamiltonian perturbations and is recovered in our theory for weak (x < 1) 
local perturbations, our analysis predicts a novel decay regime for strong (x 3> 1) 
perturbations that stems entirely from the local nature of the Hamiltonian perturbation, 
i.e. the escape-rate regime, and quantitatively describes the transition between the FGR 
and escape-rate regimes as the perturbation strength is varied. In the escape-rate regime 
the LE decays exponentially in time with a rate equal to twice the escape rate from an 
open billiard with the "hole" at the place of the perturbation. Hence the LE allows 
to mimic the decay behavior of a system without opening it. In this regime the LE 
decay rate is independent of the deformation strength x- The transition between the 
FGR regime and the escape-rate regime turns out to be non-monotonic: the rate of the 
exponential time-decay of the LE oscillates as a function of the perturbation strength. 

We would like to point out that recently there has been another study [IE] of the LE 
decay due to local perturbations. It addresses a particle moving in a two-dimensional 
array of point-like scatterers, and the perturbation of the Hamiltonian is achieved by 
slightly displacing one of the scatterers. The theoretical and experimental analysis of 
the system reveals a polynomial decay of the LE, namely M(t) ~ t~ 2 , for long times, 
which cannot be compared with our findings: that work is in the perturbative regime, 
where the eigenstates are not significatively modified by the perturbation. 

We have also performed an extensive numerical study of the LE decay to support 
our semiclassical theory. To this end we have simulated the time evolution of initially 
coherent states in the DDB. The role of the local Hamiltonian perturbation was played 
by a piston-like deformation of the billiard boundary. The results of our numerical 
simulations exhibit strong quantitative agreement with the predictions of our theory 
extended to cope with initial states given by Gaussian wave packets of a dispersion 
comparable with the de Broglie wavelength. 

While the scale separation given by Eq. ( 155]) is rather challenging to be satisfied 
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numerically it can be naturally achieved in laboratory experiments with ultra-cold atoms 
confined to optical billiards. In this work we have proposed a laboratory set-up allowing 
one to investigate the LE decay from local Hamiltonian perturbations for a wide range 
of perturbation strengths, and to observe the predicted decay regimes. We believe 
that the study of the LE decay due to local perturbations provides an example of 
physical problems for which capabilities of laboratory experiments go beyond those 
of numerical simulations. Such experiments may also reveal weak-localization-type 
quantum corrections to the LE decay which are expected from an analysis [47J of loop 
contributions [H] beyond the semiclassical diagonal approximation. 
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Appendix A. Linearization of the action integral 



Here we present the details of the expansion of the action integral (or the Hamilton 
principal function) around the central trajectory of the wave packet. This approximation 
is an important aspect of the semiclassical description of the LE (see Ref. [H] and Sec. [2] 
of the present work), and its limits require careful consideration. Thus, we devote this 
appendix to the validity condition of the linear expansion. Higher order expansions have 
been recently considered in the literature p3J EH HH] . 

As shown in Fig. [H we consider a wave packet centered at r and localized at a 
small circular region of radius a. The action integral S§(r,T',t) along a trajectory s 
starting at a point r' within this circular region at time and leading to r in a time t 
can be expanded as 



Sg(r,r',t) = S s (r,r ,t) 



+ -(r 

2 V 



ro 



dS s (v,r',t) 



d 2 S s (r,r>,t) 



(A.l) 



=r 



dr 



12 



(r' - r ) 



=r 



Here we assume that the trajectory s(r, r', t) converges to the central trajectory s(r, r , t) 
as r' approaches ro- The dot denotes the scalar (as opposed to matrix) multiplication. 
Using the identity dSg/dr' = — p§, where pg denotes the initial momentum on the 
trajectory s(r, r',t), we rewrite Eq. ( 1A.1I) as 

S s (r, r', t) = S s (r, r , t) - p s ■ (r' - r ) (A.2) 

1 , . - \®Ps , , 
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Note that in our notation — > p s as r' — > r , see Fig. [TJ In order to truncate the 
expansion (1A.2I) at the term linear in r' — r , and therefore recover Eq. ([6]), the condition 



(r' - r ) 



dps 
dr' 



=r 



< h 



(A.3) 



must be satisfied for all points r' such that |r' — ro| < o. 

To analyze Eq. flA.3j) we introduce a system of relative coordinates moving along 
the central trajectory s(r, ro,£). Thus, for any time r G [0, t] the distance between 
the phase space points (q^,p^_) and (q T ,p T ), belonging to the trajectory s(r, r',t) and 
s(r, r , t) respectively, is given by — q r = (g|, q^r) and p' T — p T = (pl,p^), where the 
superscripts || and J_ refer to the vector components parallel and perpendicular to p T . 
(Note that in the current notation q = r , q = r', q t = q[ = r, p = p s and p' = p s .) 
Then 



'dpi 




dr' _ 


r'=r \ 



dpl/dql dpl/dqfr 
dpo/dql dp^/dq^ 



(A.4) 



(<?o'> 95^=0 

For a billiard the off-diagonal partial derivatives can be neglected compared to the 
diagonal ones, so that condition flA.3|) can be replaced by 



a 



dpi dp^ 



dq\ dq^ 



(A.5) 

(9o'.?o ± )=0 

The first of the two derivatives in Eq. ( 1A.5I) is dpl/dql = —m/t for a particle of mass m 
in a billiard. To evaluate the second derivative we first linearize the trajectory s(r, r', t) 
around s(r, r ,t), so that q^r « g^g^pQ, r). Therefore, 



= dqi 



dQt 
dQo 



dqo + 



dpo 



dp, 



o ) 



(A.6) 



which leads to 



dpp 
dQo 



{dqt/dp^) c 



(A.7) 

The right hand side of Eq. flA.71) is given by the ratio of the two monodromy matrix 
elements. To facilitate our analytical presentation, we use here the monodromy matrix 
of the dynamics on Riemann surfaces of constant negative curvature 



dqildq^ dq^/dp^ 
dpi/dq^ dpi/dp^ 



cosh(At) (mA) 1 sinh(At) 
mA sinh(At) cosh(At) 



(Ai 



with A the Lyapunov exponent. For times longer than the Lyapunov time, t ^> 1/A, we 
have dpQ /dq$ ~ — mA, so that Eq. (IA.5[) can be replaced by 

a 2 mA<ft. (A.9) 
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In terms of the Lyapunov length /l — conveniently used for billiards, 

Eq. flA.91) then reads 



(A.10) 

The momentum uncertainty of a Gaussian wave packet of dispersion a is h/a, so that 
p < hj X + h/o~, with X being the rescaled de Broglie wavelength. Therefore, condition 
( 1A.10I) holds for every trajectory relevant for the wave packet propagation only if 



The action integral expansion requires condition (lA.lip to be satisfied. We finally 
note that in the limit X ^ a, which we utilize in Sec. El Eq. (lA.llj) simplifies to 
o < v^l- 

Appendix B. Piston-like boundary deformation 



w 




Figure Bl. Piston-like boundary deformation. 

In this appendix we explicitly compute as an example the length-difference function 
of Eq. fl3Tl) for a piston-like local boundary deformation, see Fig. IBll which is 
also used in our numerics. We assume that the boundary of the unperturbed billiard 
possesses a straight segment of length w that gets "lifted" by the perturbation as if an 
imaginary "piston" was pulled out. We denote the piston displacement by h. Assuming 
h much smaller than the free flight path If of the trajectory hitting the deformation, 
we treat the unperturbed and perturbed trajectories to be parallel. Then the length 
difference accumulated due to a single collision with the deformation-affected 

segment of the boundary is given by an expression analogous to Bragg's diffraction 
formula, 

« 2/icostf. (B.l) 

Here $ represents the collision angle as shown in Fig. IBll The length difference u can be 
considered independent of the collision coordinate £ for deformations such that h <w. 
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Taking into account the probability distribution function of collision angles, 
Eq. (132]) . we have (cosi9) = it/A and (cos 2 ??) =2/3. Consequently, the first two moments 
of the deformation function read 

(u) = -h and (u 2 ) = -h 2 . (B.2) 
2 3 

Similarly, the average phase factor due to a single collision of the particle with the piston 
is given by Eq. (|36|) . 
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